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THE CYLINDRICAL ANTENNA; CURRENT AND IMPEDANCE* 


BY 
RONOLD KING anp DAVID MIDDLETON 


Cruft Laboratory and the Research Laboratory of Physics, Harvard University 


1. Introduction. The definition and the determination of the impedance of a sym- 
metrical, center-driven antenna of small, circular cross section involves three major 
problems. These are first the theoretical analysis including the formulation of bound- 
ary conditions; second the apparatus and the technique of experimental measurement; 
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red 


Fig. 1. Cylindrical 
antenna with hemi- 
spherical! ends. 


and third the coordination of experiment with theory. Of these 
problems only the first is the subject of this paper; the last two 
are considered in detail elsewhere.! The present discussion is 
concerned specifically with an analytical improvement in the 
solution of the theoretical problem. 

The boundary and driving conditions in this analysis are the 
same as implied in earlier analyses,?:* and the same integral 
equation is obtained. However, the present paper introduces a 
new approach to the solution of Hallén’s integral equation in 
that it replaces a function arbitrarily chosen for reasons of 
mathematical convenience in the approximate evaluation of the 
equation by a function actually fitted to the true distribution of 
current. As a consequence, new parameters are introduced to 
replace those used by Hallén (or equally those used by Gray”) 
in the successive approximations, and as would be expected the 
resulting development shows a more rapid convergence, in so 
far as this is indicated by a relatively small difference between 
first and second order solutions. 

The antenna actually analyzed is a theoretical one in the 
sense that no exact experimental analogue can be constructed. 


Its properties are summarized as follows: 
(1) The antenna is a highly conducting cylinder of small radius @ extending un- 
broken from z= —h toz= +4 as shown in Fig. 1. Postulated inequalities are 


pa <1, a<h, (1) 


where B=q/c is the phase constant and ¢c=3 X10 m/sec. 
(2) The ends of the antenna at z= +h contribute nothing to the electrical prob- 
lem so that it is correct to write 


I,=0 at cs=Hth. (2) 


(3) The antenna is center-driven by a slice generator consisting of a disk of neg- 
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ligible thickness at the center, z=0, which is in all respects like any other piece of the 
antenna except that a scalar potential difference 


Vo = lim (p+: = $-2) = 40 — $-0 (3) 
20 


is maintained between its faces. 

(4) All other conductors and all dielectrics are sufficiently far away so that their 
individual effects are indistinguishable from the composite effect of the universe as 
a whole. If R is the distance from the center of the antenna to the nearest point on 
any other conductor or on a dielectric the following inequalities must be satisfied 


BR> 1; RD hk. (4) 


The degree in which this theoretical antenna can be realized physically is sum- 
marized briefly below. Details are found elsewhere.!45:6.7 

(a) A metal wire or rod can be constructed to satisfy completely the properties 
assumed in conjunction with (1). 

(b) If a solid cylinder with flat ends or a hollow cylinder is used (2) is not exactly 
true. A small current exists at the ends to charge the sharp edges, the end surfaces, 
or the inner surfaces of a tube near the ends. This leads to an error in / of the order 
of magnitude of a, and a consequent hidden shift in the theoretical impedance curves. 
For particular values of # near anti-resonance large errors in impedance are involved. 
A solid cylinder of length 2h along the axis with hemispherical ends as shown in Fig. 1 
is a satisfactory physical approximation that satisfies (1) and (2). 

(c) It is physically impossible to provide a slice generator. At very low frequencies 
a two-wire drive is satisfactory to approximate (3) but in this case (4) can not be 
satisfied. At high frequencies where (4) is readily satisfied a two-wire drive involves 
adjacent end surfaces and a gap in the antenna which are not taken into account in 
the theory. The effects of gap and end surfaces are compensating and may be taken 
into account roughly in comparing theoretical and experimental results by including 
a lumped capacitance in parallel with the experimentally measured impedances if the 
gap is large and a similar capacitance in parallel with the theoretical impedances if 
the gap is very narrow and the adjacent end surfaces of the antenna are very close 
together.’ A vertical antenna of length # over a conducting plane, driven from a 
coaxial line, may be a good approximation of a slice generator, but the unavailability 
of an infinite, perfectly conducting plane leads to other difficulties. 

(d) The condition for the far zone (4) can not be fulfilled at low radio frequencies 
(where accurate measurements can be made easily) because it is not possible to get 
far enough away from the earth. At high frequencies where this is possible, accurate 
measurements are difficult and the dimensions of the antenna and its driving struc- 
ture become undesirably small.} 

The analysis of the theoretical antenna subject to the conditions (1) to (4) dis- 
cussed above, can be reduced to one-dimensional form involving the total current J, 
if it is assumed that the cross-sectional distribution of the density of current is inde- 


‘L. Brillouin, Quart. Appl. Math. 1, 201 (1943). 
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pendent of the axial distribution. In effect, this means that the cross-sectional dis- 
tribution and the internal impedance per unit length may be obtained from the analy- 
sis of an infinitely long cylinder. This is an excellent approximation subject to (1). 
At high frequencies the internal impedance per unit length is given by 


1 ‘wull 
bay 49 (5) 


2ra 


with o the conductivity in mhos per meter, » the relative permeability, a the radius 
in meters, and II =42 X10-’ henry per meter. Subject to this assumption in addition 
to (1)-(4) the vector potential at any point z on the cylindrical surface of the antenna 
due to the axial current in the entire antenna is given by * 8° 

Tl h 

A,=— Izy Rye~ #812’, (6) 

an J _, 

where 


Ri = Ve — 2’)? + a, 


and J/ =J,(z’) is the axial current at 2’. 
The integral equation for the current, originally derived by Hallén,? is*® 


(7) 


4r h : 
—A,= If Ry te7#BRdz! 
Ti -h 

— j4n 


ll 


& 

gE cos Bz + 4Vo sin B| 2 | - if I(s) sin B(z — sas. (8) 
e 0 
Vo is the driving potential difference maintained by the slice generator at z=0; C; is 
a constant of integration which is later evaluated using (2); R.=cIIl =376.7 ohms 
1207 ohms. In practice the conductivity o is usually sufficiently high and therefore 
z' sufficiently small so that the last integral in (8) contributes negligibly to the final 
result." For simplicity it is omitted throughout the following analysis. If required 
it can be included readily at appropriate points with no change in the formulation. 

2. Expansion of the integral equation. In the absence of an exact solution of the 

integral equation (8) in closed form, an approximate solution may be obtained by 
expanding the integral on the left in a converging power series in terms of an ap- 
propriately chosen parameter. If a converging series is obtained and a sufficient num- 
ber of terms can be evaluated the choice of the parameter for expansion is unimpor- 
tant. If only a few terms in the series can be evaluated readily it is of great importance 
to select the parameter in such a way that convergence is so rapid that the sum of 
two or three terms gives a satisfactory approximation. The several parameters which 
have been used,®” including that introduced below, will be discussed critically and 
results compared in another paper. The general definition of all such parameters is 
formulated below. 


®R. King, Electromagnetic engineering Vol. 1, McGraw-Hill Book Co., New York, 1945, p. 241. 
9S, A. Schelkunoff, Electromagnetic waves, D. Van Nostrand Co., New York, 1943, pp. 140, 142 ff. 
0 Reference 3, equation (25). The complete derivation is given. 

4 R, King and F, G. Blake, Proc. 1.R.E. 30, 335 (1942), 

2M. C. Gray, J. Appl. Phys. 15, 61 (1944). 
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The solution of (8) may be formulated by expressing J, in terms of a convenient 
reference current such as the input current J) and a distribution function f(z) that 
is unknown. Thus let 


IT, =Inf(z); 9 Tz = Tof(z’), (9a) 
so that 
Iz = I.f(2')/f(z) = 1g, 2’). (9b) 


The relative distribution function g(z, 2’) is defined in (9b). Now let a function V(z) 
be defined by A 
V(z) = f g(z, 2’) Rete #Ridz’, (10) 
ah 
If the relative distribution function g(z, 2’) were the actual one, it would be correct 
to write J,¥(z) for the integral on the left in (8). Whatever the form of g(z, 2’), it is 
correct to write 


h h 
=A, =f IL Ryp'e- Fide’ = 1,0(z) + [yt — I,g(z, 2’)| Rye #Pidz’, (11) 
—sA —sA 

The more nearly g(z, 2’) approximates the true distribution the smaller will be the 
difference integral on the right in (11). If g(z, 2’) can be chosen accurately enough so 
that the integral on the right in (11) is considerably smaller than the term I,¥(z) for 
all values of z, it is possible to treat this term as the principal part and the difference 
integral as a correction. 

If g(z, z’) were the true relative distribution function so that the difference in- 
tegral in (11) were zero, the function Y(z) would be given by 


¥(2) 4n A, 
z= — < 
nN I, 


(12) 


That is, ¥(z) would be proportional to the ratio of the vector potential on the surface 
of the antenna at a point z divided by the total axial current at z. It is clear from (6) 
that the vector potential at a point 2-is determined largely by the current at and 
near 2, except possibly at a few points where J, is very small compared with the cur- 
rents elsewhere in the antenna. It may be assumed, therefore, that the ratio A,/I, 
is reasonably constant and predominantly real at all points along the antenna except 
at and near very small or zero values of the current. Clearly, since J, =0 at the ends 
and A, is not zero there, V(z) is infinite at z= +h. However, the product 7,¥(z) must 
remain finite and relatively small at z= +h. 

If W(z) is sensibly constant for most values of 2, it must be exactly equal to V(z0) 
at some point 2 =o, so chosen that V(zo) is a good approximation of W(z) except where 
I, is small or zero. Let 

v =| ¥(a)|, (13a) 


so that 
W(z) = Very, (13b) 


Also let a function y(z) be defined so that 
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W(z) = V+ (2), (14a) 
where 
y(z) = Vere — 1), (14b) 


If V(z) is predominantly real, (z) is a small complex correction function except at 
values of z where J, is small or zero. It is to be noted that y(z) is infinite at z= +h 
but that J,(z) is finite and small there. 

If g(z, 2’) is not the true relative distribution function but only approximate, 
(12) is also approximate, and it is still possible to write (13a, b) and (14a, b) with 
W(z) defined as in (10). Substituting (14a) in (11) and using (11) in (8) solved for J, 
in the principal term J,¥, one obtains 

[ea ths {C, cos Bz + 4Vy sin B| 2| } 
RW 


1 Ly 
_ s{tat) + [Zi — L.g(z, 2) [Reena : (15) 
—h 
This equation is exact. Like (8) it is an integral equation in the current, but the cur- 
rent appears in the integrand of a difference integral that is small. The term J,¥(z) 
is also small except near points where J, is small or vanishes, as at z= +h. 
A more useful form of (15) is obtained as follows. Let (8) be written with z=h in 


the form 
— jar 


‘ 1 i 
= {Cy cos Bh + 4V, sin Bh} - —f I} Rye 8 Finda’. (16) 
h 


RY W_ 


The term in 2‘ has been omitted in (16) just as in (15). Actually (16) is exactly equiva- 
lent to (15) when this is written with z=h. In (16) 


Ry = V(h— 2)? + a? (17) 
The desired equation is obtained by subtracting (16) from (15). It is 
— jan 
RY 


I,= {Ci[cos Bz — cos Bh] + $Vo[sin | z| — sin Bh]} 


1 h A 
- = {ito +f [72 — L.g(z, 2’) |Rrte~#Pid2’ -f I; Ryteritends . (18) 
—h aA 


This is the final exact form of the integral equation. Its principal advantage over (8) 
lies in the fact that all terms on the right involving the current are small if the rela- 
tive distribution function g(z, 2’) is correctly chosen to make the difference terms 
small. The expression (18) must be used in preference to (15) because in (18) the 
right side vanishes for all values of Bh when z= +h as required by (2), whereas the 
right side in (15) can not be made to vanish at x = +h when cos 64=0. In this case 
the arbitrary constant C, disappears from (15). 

The integral equation (18) can be expressed as the sum of a principal current 
(7,.)a consisting of the trigonometric terms and a correction current (J,), given by the 
remaining terms. The correction term (J,). can then be expanded in a power series 
in 1/¥. Thus 
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I, = (I,)o + (2). = (Tz) + Iq + (Ta) + (LT) 3 + Pe ie ce (19) 


Here (J,)., is obtained by substituting (J,)o in (72); Us), +(:)., is obtained using 
(7,)o+(2)., in (1s)e, etc. 
For convenience let 


Fn(z) — F(h) = Fue Ga(z) — Galt) = Gas (20a) 
where 
F(z) = cos Bz; Fo(k) = cos Bh; Go(z) = sin B| z|; Go(h) = sin Bh (20b) 
F,(z) = fae | ae 2) Ree BRidz! — "Pye REM Pda! — Fy1sv(z), (21a) 
-h -h 
Fy(h) = — | Fyoaer Ryle ®Pdz’. (21b) 


—h 


The first and last terms in (21a) may be combined into Fy_1,, V using (10) and (14a). 
Expressions for G,(z) and G,(k) are obtained from (21a) and (21b) by writing G for 
F throughout. 

Using (19)-(21) in (18), the complete series solution for J, may be obtained. The 
constant C; may be evaluated from (16) using (19)-(21) as described in references 
2 and 3. The resulting mth order current is 


Yo Fale) /¥"- Ys Galh)/¥* — Yo Gala)" Oo Fah) 
n= n=O nm 


j2aVo \ m0 


q, mn = 22 
(a = : (22) 
De Fah) /Y" 
n=O 
This formula may be simplified using (20a, b). The result is 
«(sin 6(2 —|2]) + D> Ma(z)/er 
J4TVo n=l 
(I2)m = ry (23) 
RY ul 
cos Bh + >> F,(h)/¥" 
where, in particular, a 
My(2) = My(z) + 7M; (2) = F,(z) sin Bh — F,(h) sin B| z| + G,(h) cos Be 
— Gi(z) cos Bh, (24) 
Ma(2) = Mi(z) + 7M: (2) = Fi(2) sin Bh — F,(h) sin B| 2 
+ G(A)Fi(z) — Gi(z)Fi(h) + Goh) cos Bz — G2(z) cos Bh, (25) 
and, as previously defined,?+.56 
PY = acme, + jas’, (26) 
With 


8 C, J. Bouwkamp, Physica 9, 609 (1942). In Bouwkamp’s paper G and Fare, respectively, the F 
and G functions in this analysis. 
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I _ TT T eg gll 
Bn = Bu + jn = Mx(0) = M,(0) + JM, (0), (27) 
the impedance of the antenna is defined by 
Zo = Vo/Io, (28) 


where Jy is given by (23) with z=0. It is 
cos Bh + > a,/¥" 


‘ —J eV nel 
(Zo)m = (Ro)m + j(X0)m = ——— (—____}.. (29) 
2n us 
sin Bh + > B,/¥* 
n=l 


This is a generalization of the formula obtained by Hallén? and others*:, as shown 
later. 

3. Functions and parameters in the Hallén solution. The expressions for the cur- 
rent (23) and for the impedance (29) depend upon the constant parameter VY, and 
this in turn depends upon the relative distribution function g(z, 2’). The definition of 
these quantities involves the following considerations: The relative distribution func- 
tion g(z, 2’) must be so chosen that it is a sufficiently good approximation of the 
actual distribution to make the difference integral in (11) small. Furthermore, it must 
be sufficiently simple in form that the integral (10) for Y(z) can be evaluated and 
separated into a principal, constant, real part | ¥(z0)| =VY and a small correction 
term y(z) as in (14a, b). 

The choice of distribution function made by Hallén depended upon the reasonable 
albeit implicit assumption that the vector potential A, at 2 depends primarily upon 
the current at and near z. If contributions from all more distant elements of current 
are small, A, may be evaluated approximately by assuming the current at all points 
to be J, and neglecting retardation. This is equivalent to setting 


gu(z, 2’) = e78R1, (30) 


The subscript H will be used to designate parameters and functions in the Hallén 
analysis. With (30), (10) gives 


* da! . h+z . h—32 
Vu(z) = — = sinh" + sinh™! : (31) 
—h Ry a a 
Alternatively and equivalently 
2? 
Vi(z) = Q2+41n (: _ 5) + 8(z), (32) 
where 
2h 
Q = Vx(z = 0) = 2iIn—,; : (33) 
a 


om Weg WG} « 


If the function Wy(z) in (31) is plotted as a function of z/a for a range of values of 
the ratio k/a, it is found to be moderately constant for large ratios h/a except with 
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z near h. Subject to the condition a?<h?, the maximum value of Vy(z) is Q, the aver- 
age value is Q—2+42 In 2=0Q—0.614, the value at z= +h is }0+1n 2. For Q>15 the 
average value or the maximum value are satisfactory approximations. In view of the 
fact that Vx(z) becomes smaller at z= +h instead of becoming infinite as it would if 
the correct distribution function were used, it is clear that the curvature of Vx(z) 
is the reverse of what it should be. Therefore, the maximum value Q is probably the 
best approximation of Vy(z) and this was Hallén’s, although not explicitly for this 
reason. Thus the Hallén analysis sets 


2h 
W(Zo) =vy =2=2ln—, (35a) 
a 


vate ein (: a = + a(z). (35b) 


The Hallén expressions for the current and the impedance are given by (23) and (29) 
with Q written for Y and with appropriately modified functions F,(z), Fa(h), Ga(z), 
and G,(h), n>0. The functions with »=0 are independent of the choice of VY. The 
Hallén functions are 


h e7 BR, 
Putin Seee f ae — ae (36a) 
= Ry 
h e7 Rin 
Fasthye Ss f Ph de. (36b) 
~h Ry 


Gan(z) and Gan(h) are obtained from the above by writing G for F. These functions 
have been evaluated elsewhere®’ for n=1 and »=2. The first order distribution of 
current and the first order impedance have been calculated and represented graphi- 
cally*:1!; the second order impedance has been evaluated by Bouwkamp.™ The Hal- 
lén formula for the mth order current is 


suas sin 6(h —|2|) + 2 Mau(z)/Q* 
(12) mH = 7 : = ’ (37) 


RQ m 
cos Bh + >> Fan(h)/Q* 


n=1 


cos Bh + > any/2Q" 
c n=1 


Ou ~ (38a) 
sin Bh + >> Ban/Q" 


n=l 


Here 
I Ir I II 
AnH = Ayn + Jan = Fin h), (38b) Ban = B, + IB. = Man(0). (38c) 
The functions a; and §; are tabulated and represented graphically in references 2, 3, 
11; the functions a and #2 as calculated by Bouwkamp* using graphical methods are 
listed in Table I and plotted in Figs. 2 and 3. 
4. Functions and parameters in the improved solution. The relative distribution 
function g(z, 2’) in (30) is the simplest and the most obvious one if an attempt is 
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TABLE I 

- Bh al ad! 6 pi 

0 0 0 0 0 

0.2 —0.16 0.03 3.07 = 

0.4 —0.53 0.13 5.20 0.03 
0.6 —1.07 0.39 6.50 0.24 
0.8 —1.67 0.80 7.14 0.78 
1.0 —2.17 1.31 6.78 1.74 
1.2 —2.66 1.84 5.48 3.04 
1.4 —3.00 2.31 3.34 4.97 
1.6 —3.23 2.73 0.45 7.06 
1.8 —3.34 3.04 —3.06 9.33 
2.0 —3.33 3.30 -7.03 11.81 
2.2 —3.16 3.48 —11.25 13.98 
2.4 | —2.84 3.58 —16.22 16.08 
2.6 —2.35 3.58 —20.83 17.28 
2.8 —1.59 3.40 —24.71 17.72 
3.0 —0.61 2.99 —27.54 17.50 
3.2 0.50 2.27 —29.02 16.85 
3.4 1.58 1.37 —29.29 15.75 
3.6 2.59 0.28 —28.30 13 .84 
3.8 3.49 —0.83 —26.35 11.34 
4.0 4.33 —2.00 —23 .38 8.28 
4.2 5.03 —3.09 —19.60 4.73 
4.4 5.41 —4.13 —15.14 0.21 
4.6 5.46 —5.04 —10.26 —4.84 
4.8 5.20 —5.67 —4.21 —10.09 
5.0 4.66 —6.08 +2.41 —15.10 


made to solve the original integral equation (8) as was done by Hallén. On the other 
hand, if the formal solution is carried through to obtain (23) without previously 
selecting g(z, 2’) as has been done in the present analysis, it is perfectly clear that the 
leading term in the distribution of current for any value of the distribution function 
must be of the form 

I, = Kfi(z) (39a) 
with 

filz) = sin B(k — | 2] ). (39b) 

K is an amplitude factor independent of z. Accordingly, an approximate relative dis- 
tribution function is 

sinB(h—|2'|) fale’) 


sin B(h — | 2] ) 7 fi(z) 


This function is known to be a very much better approximation of the actual current 
then the function assumed by Hallén, gu(z, 2’) =e#*. The function fi(z) =sin B(A — || ) 
actually is proportional to the principal part of the current; the function e%®® 15 not. 
Using (40) and (10) we obtain 


gri(z, 2’) > (40) 
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Fic, 3. The parameters 8 and gi as a function of fh. 
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A 
Wri(z) = eri(z, 2) Role BR iz’, (41) 
—h 
This function involves the factor f,(z) =sin B(h—|z|) in the denominator of the 
integrand. Since this is not a function of z’ it is a constant in the integration. There- 
fore, it is convenient to introduce the function 


¥ilz) = f * f(a!) Rete Bada! (42a) 
so that ci ia} 
1(z 
Wrilz) = fle . (42b) 
The function y;(z) can be written in the form 
Vi(z) = C(z) sin Bh — S(z) cos Bh, (43) 
where 
C(z) = f oe Ba’ Rte #BRigz’, (44) 
‘ 
S(z) = f au B| 2! | Rye #Rigz’, (45) 


These integrals are evaluated in the Appendix both in general and in a simpler ap- 
proximate form. The latter is a good approximation if, as assumed throughout this 
analysis, h*>>a*. Curves for C(z) and S(z) as calculated using the simpler forms which 
apply in this analysis are given in Figs. 4-7, 20-23 for Bh=a/2 and zw and for Q 
=2 In (2h/a)=10 and 20. It is to be noted that 


Wiz) =C(z); Bh = «/2, (46) 
¥ilz) = S(z); Bh =m. (47) 


It follows that the plots of C(z) with 64 =7/2 are also plots of ¥:(z); these are given 
in Figs. 4 and 5 for 2=10 and 20. Similarly plots of S(z) with Bh=7 are also plots 
of ¥,(z); these are given in Figs. 6 and 7 for 2=10 and 20. The function y(z) is seen 
to have a very small imaginary part so that it and Vx,(z) =¥;(z)/fi(z) are predomi- 
nantly real, in confirmation of the assumption made in conjunction with (14). Ac- 
cordingly, the parameter YW = | W(z0)| ‘defined in (13a) may be chosen to be 


Wri =| ¥x1(0)| =|¥(0)|; Bh = 2/2; (48) 
Wri =| Valk —2/4)| =[vilh-2/4)|; Bh =o. (49) 
The function ie | 
1\% 
| Wxi(z)| = rye (50) 


is plotted in Figs. 4-7. For Bk =1/2 and both for Q=10 and 20 it is seen to be quite 
constant over the entire length of the antenna except near the ends where it becomes 
infinite, as it should. For Bh=7 the function becomes infinite not only at the ends 
but also at the center. The infinity at the center is a result of approximating the 
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current by the distribution function f,(z) =sin B(h—|2| ). With Bh=a7 and z=0, this 
vanishes so that Wxi(z) necessarily diverges. Unlike the infinity at the ends, the in- 
finity at z=0 is due to the fact that fi(z) and hence gx(z, 2’) are approximate and not 
exact distribution functions. Actually, the current does not vanish at z=0; it merely is 
small so that W(z) does not become infinite. The fact that J, is small at and near 
z=0 does not mean that ¥(z) necessarily becomes very large. ¥(z) is by definition 
proportional to the ratio A,/Z,, and A, is determined largely by the current at z. 
Hence A,mo is small if Z,.9 is, and the ratio may remain moderately constant. Fur- 
theremore, since A, at z=4—}/4 is determined principally by the large (near maxi- 
mum) currents at and near z=h4—)/4, it is affected only very slightly by a small 
current at z=0. Therefore A, at z=h—d/4 and Vxi(h—2/4) will not be sensibly dif- 
ferent if a fictitious zero current is assumed at z=0 or.an actual small current. Ac- 
cordingly the function |Wailh —/4| is a good approximation of Vxi(z) for the actual 
current everywhere (including z=0) except near the ends, z= +h. 

Although the qualitative argument to show that Vxi(z) is sensibly constant and 
approximately equal to | Wxi(h—2/4)| for all values of z except the ends is sound, 
it can be verified directly using Hallén’s first order distribution. It has been shown? 
that a very satisfactory approximation of the Hallén first order current is given by 


(Se) ane sin p(h—[s|); ~s 8h <2m, (51) 
1 — cos Bk 2 


I, = 13’ 
where If’ is the component of current at z=0 in phase with the driving potential 
difference and I, is the maximum value of the component of current in phase quad- 
rature with the driving potential difference. Jn occurs at z=h—d/4. With 


k=Ié'/In; |r| <4, (52) 
it is possible to write (51) in the form 
: ; ., {cos Bz — cos Bh . 
I, = jIm {sin B(h — | z|) - a(S) = jl m f2(2). (53) 
1 — cos Bh 


With this approximate current, an appropriate distribution function g(z, 2’) is defined 
by 


I! sin B(h — | 2! | ) — jR(cos Bz’ — cos Bk)/(1 — cos Bh) 


(2.5) =< —— = : 
I, sin B(k — | 2| ) — jk(cos Bz — cos Bh)/(1 — cos Bh) 
_ fle) oe 
=o : S Bh < 2m. (54) 


The ratio factor & is negative and small compared with unity. It is plotted in Fig. 8 

as a function of 8h from the data of Figs. 9-11 in reference 3. Only values of BA near 

w are used because for 8A not near integral multiples of a the distribution (cos Bz 

—cos Bh) does not differ greatly from sin B(h — | 2| ). At 8h =2/2 they are identical. 
Using the notation (42a, b), 


A 
¥2(2) = J fo(z’)Rrve PF 1dz’, (SSa) 
—h 2 
so that 
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¥2(z) ; 


Wx2(z) = (SS5b) 
. g fo(z) 
The function yW2(z) can be written 
Ya(s) = C(s) sin Bh — S(z) cos Bh — jk [See (56) 
1 — cos Bh 
where C(z) and S(z) are defined in (44) and (45), and E(z) is given by 
A 

E(z) = Ry e— PRidz’, (57) 


—h 


Fic. 8. The quantity k=I$’ /In asa function of 8h near anti-resonance. 


This function is evaluated in the Appendix both in general and in a simpler ap- 
proximate form valid when it is possible to write a?<h?, as in the present analysis. 
E(z) is plotted in Figs. 24 and 25 for Bh=7/2 and w and with 2=10 and 20. The 
function y2(z) is necessarily predominantly real because it is known that the first 
two terms in (56)—these are identically ¥,(z)—are predominantly real and that k is 
small. The functions | ¥2(z)| and | Vxe(z)| = | ¥2(z)| /| fo(z)| are shown in Figs. 9 and 
10 for Bh=am and 2=10 and 20. It is seen that | Wx2(z)| does not become infinite at 
z=0, and is reasonably constant and equal to |Wx2(h—)/4)| = | Wa(h —)/4)| for all 
values of z except at the ends where it becomes infinite, as it should. Comparison of 
Figs. 9 and 10 with 6 and 7 shown that | Wo(h —)/4)| differs only slightly from 
vi(h—/4)|. The difference is greater for the smaller value of Q. It follows that 
¥i(h—X/4)| is a satisfactory parameter even for Bh=7. If desired | W2(h—2/4)| may 
be used especially for small values of 2, but the difference is not over about 3% for 
2210. 
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The conclusions of the above analysis may be generalized and summarized as 

follows: 1. The relative distribution function 

sin B(h — | 2’| ) 

sin B(h — ||) 


is a good approximation for all values of #. 2. Suitable parameters for expansion are 


gix(2, 2!) = 


[f= 20s 


Fic. 11. The parameter y, Eqs. (48) and (49), as a function of Bh, 
for 2=10, 15, 20. 


|Wx1(0)| =|y1(0)| for Bh Sx/2; |Wxi(h—D/4)| =|ya(kh-V/4)| for Bh2 2/2. The fol- 
lowing notation will be used from here on 

| vxx(0) | = | ¥1(0) |; Bh S 1/2 
| ¥xi(h — 0/4)| =| vi(h -— 2/4) |; BAS x/2. 
Since Vxi(z) has such a small imaginary part and is so well represented by (58) ex- 
cept at the ends of the antenna, the correction function y(z) in (14a) is sufficiently 


small to be neglected. 
The parameter y as defined in (58) is plotted as a function of Bh for 2=10, 15, 20 


yaw-| (58) 
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in Fig. 11. For Q=10 the curve in solid line for Bh> 7/2 is PAC —\/4)|, the curve 
in broken line is | wilh —~d/4)]. The curves for 2=15, and 20 with Bh>7/2 are 
|ya(h—2/4)]. 

5. Distribution of current. The distribution of current is given by (23), the im- 
pedance by (29) with ¥(=y) obtained from (58) or Fig. 11 for the value of A/a in 
question. The functions Fo(z), Fo(#), Go(z), and Go(k) are unchanged; they are de- 
fined by (20b). Upon substituting (40) in (10) and using (10) in the form (14a) in 
(21a) with »=1, this becomes 


h 
Fix(z) = Fo x1(z) ocd Fo, Ry 'e7 ®F1dz’, (59) 


ah 


Upon comparing (59) with (36a) written for »=1, it follows that 


Fix(2) = Fin(z) + (¥ — 2)Pos) x. (60) 
Since ¥(z) is not involved in (21b), 
Fix(h) = Fin(h). (61) 


Upon substituting (60) in (21a) with y(z) =0 and »=2, this becomes 


h 
Fox(z) = (Fis)xWxi(z) — | (Fie) x Rr'e~#*1dz’. (62) 
—A 
Using (60) and (61) in (62), this gives 
4 4 
Fox(z) = (Fix — (Fie) x Rrtei#R ids’ — (y — Q) Fos Rye ®Fidz’, (63) 
—h ~h 


Upon comparing (63) with (36a), this time written with 2=2, and using (60) and 
(61) as well as (36a) written with »=1, we see that 


Fox(z) = (Fis)a¥ + VY — 2) (Pos) + Fan(z) ~ (Fic) xO 
| — (¥ — 2)[Fos)a¥ — Fin(z) — (Y — 2)(Foe) x]. 
When terms are collected there results, 
Fox(z) = Fon(z) + (W — Q) (Fiske + WY — 2DFin(z) + ( — 2)2(Fos) x. (64) 
Using (21b) with (60) and (61), we find 


Fox(h) = Fan(h) + (Y — Q)Fin(A). (65) 
Subtracting (65) from (64), we have 
(Fa:)x = (For)x + 2h — Q)(Fic)n + (Y — 2)?(Fos) x. (66) 
Repetition of the above procedure to evaluate (F,,.)x leads to: 
(Fandx = 3 = on W-— DF W)n; =n 20. (67) 


Expressions for Gex(z), Gex(h), (Gez)x, and (Grs)x are obtained from (64)-(67) by 
writing G for F. 
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If the functions F and G are combined to form the functions M as defined in (24) 
and (25) for the first two terms, the result is 


Mix(z) = Min(z) + (¥ — Q) sin (A — [2] ), (68) 
Mox(z) = Mon(z) + 2(W — 2)Min(z) + (Wy — b)? sin B(A — | 2] ). (69) 
In general, 
i n! 
Ma x(z) = 2) ——— (v — b)*Mnin(z); 12 2 0, (70) 
imo (n — i)!1! 
where it is understood that 
Mon(z) = sin B(k —| 2|). (71) 
Similarly, 
azl (n — 1)! ; 
Fax(h) = >, (vy — Q)Fy-in(A); n= 1. (72) 


wo (2 —i— l)id! 


Upon substituting (70)-(72) in (23) the general expression for the mth order 
current becomes 


(Dam sin Bh —|2|) + 5 Dasa) Maw (2)/0" 
janVo n=l 


(Iz) m = (73) 
Ry * 
cos Bh + >> (Dn) m—Fan(h)/y" 
where with oe . 
x= 5 (74) 
the D’s have the following significance: 
IS a en eR ae eh 
(Di)m = 1 + 2% + 2? + x8 ape at Be 
(D2) m = 1 + 2x + 3x? + 4x8 +:- 
(Ds) m = 1 + 3x + 6x? | +-- (75) 
(Dim = 1 +42 | +: 
(Ds) m = 1 oe 


It is interesting and significant to note than when the series in (75) are summed 
for an infinite number of terms, i.e., m—> ©, then 


qr} 1 1 v\" 
adx*?\1— x (1 — x)" 2 


With these values of (D,)m and an infinite number of terms, (73) is identical with the 
expression (37) obtained by Hallén. Furthermore, if a—0, Y—@ for all values of 8h, 
so that (73) approaches (37) as the radius a approaches zero. 

It is important to note that if a finite number of terms is used in (73), all terms 
belonging to a given order m of solution must be retained and no others. That is, if 
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Fic. 13. The quantity D, for the first order theory as a function 
of Bh, 2=10, 15, 20. 


Fic. 12. The quantities D, and D, for the second order theory as a 
function of Bh, 2=10, 15, 20. 
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KING~ MIDDLE TON 


04 0. 08 Lo 


Fic. 14. First order current for 8h=2/2 in units of V3/6002Dx. 


KING- MIDDLE TON 
HALLEN 


“140 “1.20 “1.00 “08 E . ° 0.2 04 


Fic, 15. First order current for 6h =3x/4 in units of Vj/600Dz. 
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an mth order solution is evaluated, only terms contributed by Myx(z) and Fyx(h) 


with »=0, 1, 2,--+, m are used. It is readily verified that this is equivalent to 
writing 
(D,) 1+(1 ~)+(1 ay 
1/2 y y ’ 
Q 
(D2)2=1+ 2(1 = =) (77a) 
(Ds)2 = 1. 


Fig. 16. First order current for Bh=- in units of Vj/602Dz. 


The expressions (D1)2 and (Dz): are plotted as functions of Bh for 2=10, 15, 20, in 
Fig. 12. Similarly it is correct to set 


(D,)o = 1, 

Q 
(D:)1 =1+ (1 - =) (77b) 
(D2), = 1. 


The function (D,); is shown in Fig. 13 for Q=10, 15, 20. 

The first order distribution of current as calculated on the one hand from the 
Hallén formula (38) in reference 3, and on the other hand, from (73), are shown in 
Figs. 14-16 for Bh=7/2, 3x/4, w, and with 2=10, 20. The function f” and f’ when 
(73) is written in the form 


nV 
R.OQDy 


(I,)1 = (f" + Gf) (78) 
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are plotted in the figures. Numerical values of Dy are given in reference 3 where D 
is written instead of Dy. 

Apart from the change in the input current which is discussed below in terms of 
the impedance, the general shape of the two sets of curves is much the same. The 
new, more exact theory leads to a distribution with somewhat greater relative ampli- 
tudes nearer the outer parts of the antenna, and with a somewhat larger component 
in phase with the driving potential difference. For Bh = the first order solution of the 
new theory is the same as the first order solution of Hallén’s theory if y is written for 
Q. Since y is less than Q, this means the new first order distribution is the same as 
Hallén’s first order distribution for an antenna of greater radius, but only for Bh=7. 

6. The impedance. The formula for the impedance according to the new, more 
exact theory is 


F cos Bh + >> (Dn) m—10n/P” 
= e n=} 


Ves 2a 


a $< (79) 
(D1) m sin Bh + > (Dn41) mBa/¥* 


n=l 


where a, and §, are defined in (38b, c); a; is tabulated in reference 11; a2 in Table I. 
Curves for (Ro)m and (Xo)m as calculated from (79) are given in Figs. 17-19. Both 
second and first order solutions are shown for 2=10, 15, 20. These are calculated 
from 


A, + jA2 


Zo)m = 6 
(Zo) WB + iB: 


ef (tan1A2/A1—tan1B2/B3) (80) 


where for the second order solution 


A; = (D;) an /p + (D2) 12 /¥', 
Az = — [cos Bh + (D:):01/¥ + (D2)s02/¥'], 
By = (Ds)2 sin Bh + (D2) o8i/¥ + B2/y", 


BE OS. (81) 


with (D,)2 and (D2). given by (77a) and (D,);, (D2): by (77b). 
For the first order solution 


A, = (Dy)eas 5 


A, = — [¥ cos Bh — (Doar), 
By = ¥(Di), sin Bh + BL ; 
By = Bi, (82) 


where (D1). and (Dj): are given by (77b). 

The impedance calculated according to the new, more exact theory differ con- 
siderably in some details but not in major outline from that obtained from the 
Hallén theory as calculated by King and Blake,'! King and Harrison,’ and Bouw- 
kamp.™ In general, resistances at antiresonance are smaller and occur at smaller 
values of 8h; resistances at resonances are greater and likewise occur at smaller 
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values of Bh. These differences are most significant for large values of the radius of 
the antenna. A critical discussion of impedance calculated from the theory here 
presented and comparison with experiment and with the theories of Hallén,?* Gray, 
Schelkunoff,*:' and others is reserved for a sequel to this paper. Accuracy of the re- 
sults and convergence of the series involved also will be discussed therein. 


Fic. 19. The input resistance and reactance of a very thin, cylindrical, 
center-driven antenna, 2=20, h/a=1.1 X10, for the first and second order 
theories. 


Using a specially constructed coaxial line and driving conditions that approxi- 
mate as closely as possible the idealized slice generator D. D. King has measured 
the impedance of cylindrical antennas with hemispherical ends. A complete descrip- 
tion of the apparatus, of the method, and of the results will be contained in a doctoral 
dissertation and in a paper to be published in another journal. A cross-section of the 
results involving all of the critical values for 2=10 are shown below together with 
the corresponding theoretical results of the theory outlined above. The agreement is 
seen to be good for a/J quantities in the second order theory, only approximate in the 
first order theory. 
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TABLE II 
Ants [X [ mia Resonant ball bal if 
resonant « —Bhanti-res. nen — Bhree. |Roat Bh=—|Xoat Bh=— 
R, |X| max Ro 2 2 2 
Experimental 
Results by 800 1.95 .60 71.5 .098 85 47 
D. D. King e 
Theoretical 
Results 860 1.80 61 71.0 .094 88 42.5 
Second Order 
Theoretical 
Results 840 1.27 .39 64.8 .065 67 30 
First Order 


APPENDIX: INTEGRAL FUNCTIONS 
The Functions C(z) and S(z). 


h 


C(z) = Rite“ ?*1 cos Ba'dz! =f (Ryte— 81 + Ry le- HRs) cos Bz'dz’, (1) 


aA 


S(z) =f Ry'e—#8* sin B| 2! | dz’ = f (Ryte7 Ft + Ryte— Rs) sin Bz’dz’, (2) 
~A ‘ C1) 


where 
= V/(z — 2’)? + a? (3a) Rz = V(z + 2’)? + a2. (3b) 
These integrals can be written in the form °- 
Ci) = $i ti t+ls+ 1), (4) 
— J : 
s@) = 2 n-ne n-mW=Fh-nth-h) (5) 
where. . 
Ih= f Ryle BRi-2 de! = eibe f Ryle BRrte-2) gy’ (6) 
0 
h 
Ih = f Rp len iB (Rite ide! = en ibe f Ryten iB Rist 202! (7) 
h 
Ts =f Rp e7 iB(Ra-2 de! = e- mf Rze- BR a2) dz! (8) 
= f Ryle—BRrt#) do! = eibs f Ryle" BBrtete ds! (9) 


The four integrals (6)—(9) can all be reduced to the form 
% 
f vle-*do (10) 
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by making the changes in the variable and in the upper and lower limits listed in 


Table III. 
: TABLE II] 
dy 
Integral w vo for 2’=0 | v forz’=h 
. r PT PT eT , ° aims J . ° 
hh | je(Rt2~2) =$elv@ ara (e'—2)]] ja[ 1] = | jatets) | 9am) 
oy 1 
. . re ee a ae *, a—s v ° ° 
I; jB(Ri—2-+2') =j8[V/ (2 —2)?F0?+(2’—2) ]| 78 [en] =F. G8(Ro—2) | j8(Rutun) 
. eos 7 «pf 2 +2 ] ; F 
Is j6(Ra—2—2') =j8[./(2’+2)*-Fat—(2’+2)]| 78 [1] =— 5 G8(Ro—z) | j8(Ra—t) 
2 2 
: { ., [2 +2 9 : . 
Te | j8(Retets) =elV eT C'+9)| jo[ 41] = | jptRets) | §8(Rutu) 
2 2 
Ro=V2*+03; t=ht+e2; m=h—z; Ryu=+/u2+a7; Ru=Vui+a?* 


In terms of exponential, sine, and cosine integrals, 


f qd e7" 
ja 


du = Ei (— jb) — Ei (— ja) = Ci (b) — Ci (a) —j Si(b) +7 Si(@). (11) 


u 


With (11) the several integrals (6)—(9) may be expressed as follows in terms of the 

exponential integral and the sine and cosine integrals 

I, = — e#*{ Ei [— j8(Ris — «1)] — Eir[— j8(Ro + 2)]} (12a) 
= — e+! CiB(Ria — 1) — CiB(Ro + 2) — f SiB(Ris — u) +7 SiB(Ro + 2)}. (12b) 


I, = Ce 


Ts 


e~ 


e~ 


{Ei [— j8(Ria + u:)] — Ei [— j8(Ro — 2)]} (13a) 
82! CiB(Ris + ui) — CiB(Ro — 2) — jf SiB(Rin + ws) + 7SiB(Ro— 2)}. (13b) 
62{ Ei [— j8(Rea — u2)] — Ei [— j8(Ro — 2)]} (14a) 


= e~##{ Ci B(Rex — u2) — CiB(Ro — 2) — 7 SiB(Ren — ue) + fSiB(Ro — 2)}. (14b) 


I 


i 


e+! Ei [— j8(Rea + u2)] — Ei [— j8(Ro + 2)]} (15a) 
e#*{ Ci B(Rea + us) — Ci B(Ro + 2) — 7 SiB(Rea + m2) +7 SiB(Ro+2)}. (15d) 


Upon combining the several integrals according to (4) and (5), 


C(z) = 


C(z) = 


ge%*{ Ei [— j8(Res + u2)] — Ei [— j8(Rus — «)]} 
+ $e-#*{ Ei [— 78(Ris + m)] — Ei [— j8(Rea — u2)]}. (16a) 
dei8e{ CiB( Re, + ue) — CiB(Ris — m1) — 7 SiB( Rea + m2) + 7 SiB(Ris — u1)} 
+ $e-#2{ CiB(Rix + m1) — Ci B(Rea — ue) — 7 SiB(Ria + 1) 

+7 SiB(Re — us)}. (16d) 
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S(z) = = ei {Bi [— j8(Res + u2)] + Ei [— j8(Ris — m1)] — 2 Ei [— j8(Ro + 2)]} 
Ea SE [i8(Ris + u1)] + Ei [—J8(Rea — u2)] — 2 Ei [—78(Ro — 2)] }. (17a) 


S(2) = = ei[Ci B(Ren + 2) + Ci B(Ris — 1) — 7 SiB(Ron + 142) 
— 7 SiB(Ris — m:) — 2 Ci B(Ro + 2) + j2 SiB(Ro + 2)] 
ft = e~#[Ci B(Rin + 1) + Ci B(Ren — 2) — 7 SiB(Ris + m1) 


— j SiB(Rex — m2) — 2 Ci B(Ro — 2) + j2 Si B(Ro — 2)]. (17b) 
In trigonometric form 
C(z) = $ cos Bz[Ci B(Rea + m2) + Ci B(Ria + 1) — Ci B(Rex — m2) 
— Cip(Ris — m1) — 7 SiB( Rea + ue) — 7 SiB(Rin + m1) + 7 SiB(Res — ue) 
+ 7 Si B(Ru — m1) 
+ 4 sin Bz[Si B(Rea + u2) — SiB(Rix + m1) + Si B(Ren — u2) 
— Si (Ris — #1) + 7 Ci B(Rea + ue) — j CiB(Ris + m1) +7 Ci B(Rea — m2) 
—jCiB(Ri,—m)]. (18) 
S(z) = 4 cos Bz[Si B(Rex + m2) + Si B(Ris + m1) + Si B(Ren — u2) 
+ Si B(Ri, — u1) — 2 Si B(Ro + 2) — 2 Si B(Ro — z) + 7 Ci (Ren + ue) 
+ 7 CiB(Rin + m1) + 7 Ci B(Rea + ue) + 7 Ci B(Rin — m1) 
— j2 Ci B(Ro + 2) — j2 Ci B(Ro — 2)] 
— 3 sin Bz[Ci B(Res + u2) — Ci B(Rin + m1) — Ci B(Rex — 2) 
+ Ci B(Ri, — 41) — 2 Ci B(Ro + 2) + 2 Ci B(Ro — 2) 
— J SiB(Rea + ue) +7 Si B(Ris + ui) +7 Si B(Rea — ue) 
— j SiB(Ris — m1) + 2 Si B(Ro + 2) — 72 Si B(Ro — 2)]. (19) 


R=V2P2+0; w=ht+s; m=h—2; Ryn = V2 + a; Ruy = Vie + a? 


These are exact expressions for the integrals (1) and (2). They are valid for all values 
of the argument z and the parameters h/ and a. 


The integral E(z). ‘ ia 
e? 1 
El) = f dat. (20) 
- 


Let the variable be changed by setting 
BR, = BV (2 — 2’)? + a? = VU? + V?, (21) 
U = Bz — 2’); V = Ba. (22) 


where 
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The integral then becomes 


BR» 
E(z) -f (U2 + V2)-M2e- FV )1/2GU 
p 


Ry 


Rep 
= f (U? + V2)-1/2 cos (U? + V%)1/2dU 
8. 


Rin 
BR 
a if (U2? + V2)-12 sin JU? + V2)12dU (23) 
aR 
with . 
Ro = V(h + 2)? + a; Ru = V(h = 2)? + a. (24) 
Let the following symbols be introduced: 
Cuv x =f (U? + V*)-1/2 cos (U2? + V%)1/2dU, (25) 
0 
Suv x -f (U? + V2)-1/2 sin (U? + V2)!/4dU, (26) 
0 
.These functions satisfy the conditions Cuv (—x)=—Cuv x, Suv (—x) =—Suv «x. 
In terms of the notation (25) and (26), the integral (20) becomes 
E(z) = Cuv BRe, — Cuv BRis — 7 Suv BRea + 7 Suv BRis. (27) 


This is an exact expression for the integral (20). 
Approximate Expressions for C(z), S(z), and E(z). 


If the parameter a appearing in Ra=W/(h+2)?+a?, Ru=~V(h—z)?+a? is small 
compared with h, useful approximate expressions for the integrals C(z), S(z), and 
E(z) may be derived as follows. Expanding the integral cosine using 


Ciz=C+inz—Ciz= fo w(t — cos w)du (28) 
0 
where C is Euler’s constant, one obtains 
Ci B(Rex — u2) = C + In (Rex — uz) — Ci B(Rea — m2), (29) 
Ci B(Rin — #1) = C + In B(Ris — 1) — Ci B(Rin — m1), (30) 
Ci B(Ro — z) = C + In B(Ro — 2) — Ci B(Ro — 2). (31) 
However, 
Ray — 2 = V(h + 2)? + a? — (hk +2), (32) 
Ru — uy = V(h — 2)? + a? — (h- 2), (33) 
Ro- 2 = V2? +a? — 2, (34) 


are all of order of magnitude a, so that the arguments B(Ro,—u2) and B(Ru—w) are 
of magnitude 8a. It follows that since the arguments of Ci B(Re.— 2), Ci B(Ris— 1) 
and Ci 8(Ro—z) are small these functions may be expanded in series. The leading 
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term in each case is of the order of magnitude 38a? so that the Ci terms in (29)-(31) 
are negligible compared with the logarithm. Hence, 


Ci B( Ren = Ue) =C + In B(Ror = U2), (35) 
Ci B(Ris — #1) = C + In B(Ria — uy), (36) 
Ci B(Ro — z) = C + In B(Ro — 2). (37) 


Since the functions Si B(Ru~—) and Si B(Ra,—wu2) are of order of magnitude 8a, 
they are negligible compared with Si B(Re,+wue) and Si B(Ri,+41) except very near 
the ends z= +h. If Ci B(Rea+u:2), Ci (Ritu), and Ci B(Ro+z) are expanded using 
(28) and the relation 


2 2\1/2 
ee ee (38) 
a a 
and the approximations, 
CiB(Rea + 2) + Ci 2Bu2 = Ci 28(4 + 2), (39) 
CiB(Ria + ur) = Ci 28a, = Ci 26(4 — 2), (40) 
CiB(Ro + 2) = Ci 262, (41) 


(18) and (19) reduce to the following approximate forms: 


C(z) = — }cos Bz [Ci 28(h + 2) + Ci2B(A — 2) + 7 Si2B(A + 2) + 7 Si28(h — 2)] 
+ }sin Bz[Si28(h + 2) — Si2B(h — 2) — 7 Ci2B(k +2) +7 Ci26(h — 2)] 


h h- 
+ cos Bz [sinn- 2 . + sinh- ‘|. (42) 
a a 


S(z) = 4 cos Bz[Si 28(h + 2) + Si 28(h — z) — 2 Si 28] 2} 
— 7 Ci 2B(h + 2) — 7 Ci 2B(hk — 2) + 27 Ci 282] 
+ } sin Bz [Ci 28(h + 2) — Ci 26(h — 2) + 7 Si 28(h + 2) 
— 7 Si 2B(h — 2) — 27 Si 26z] — sin B| z| Ci 262 


+ sin p| z| [ sik eS a ginhee "| 
a 
— 2sin Bl | [sist ai — sinh7! ily (43) 
a 


The last factor in (43) is written in the expanded form shown in order to contain 
Wy(z) =sinh—1(h+2/a)+sinh—1(h—z/a). The remaining two terms may be written in 
the following approximate form if desired, 

ning WLI ae h+|elt [(e+2)? + a2)!” 
— | sinh7? ———— — sjnh-! —— | = — In 
[ 2] + (@? + a%)1? 


a a 


= In (=) (44) 


The function (42) is shown graphically in Figs. 20 and 21 for Bh=7; in Figs. 4 and 
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5 for Bh=2/2. The function (43) is in Figs. 22 and 23 for Bh=7/2; in Figs. 6 and 
7 for Bh=7. 

An approximate expression for E(z) is obtained from (23) by adding and sub- 
tracting 


BR2b 
f (U2 + V2)-V7aU, 
8B. 


Rin 


8 

+ (U2? + V2) dU, (45) 

BRiy 

If the small quantity V=8a is neglected in the first two integrals in (45) these re- 
main everywhere finite and vanish at U=0. This is not true in the last integral in 
which Y plays an important part. If V is neglected in the first two integrals both in 
the integrand and in the limits, but retained in the last integral, the following ap- 
proximate expression is obtained: 


B(eth) B(eth) 
E(z) = -f | U|-(1 — cos U)dU — if U-! sin UdU 
a 8 


(s~h) (s—-h) 


BRab 
+ (U2 + V)"4dU. (46) 
BRin 
Because the magnitude of U and not U itself appears in the denominator of the first 
integral, this must be evaluated in two steps for the ranges 2’>z and 2’ <z. It is not 
necessary to write | U| in the second integral because the integrand does not change 
sign as 2’ passes through z. Hence with U =8(z’—z) = — U the first integral in (46) 


becomes 
Q 


B(eth) = ae 
-f U-(1 — cos U)dU — U-\(1 — cos U)dU 
0 B(s—h) 


B(h+2) 
= -f U-(1 — cos U)dU 
0 


B(A—z) = ees gh 
-f U-"(1 — cos U)dU. (47) 
0 


Using (28) in (47), the sine integral in the second integral in (46), and evaluating the 
third integral directly, we write (46) in the form: 


E(z) = — Ci@(h + 2) — Ci B(h — 2) — j SiB(h + 2) — j SiB(h — 2) 


h = 
+ sink ( ei *) + eink! ( ‘). (48) 
a a 


Use has been made of the fact that Si 8(z—h) = —Si B(h—2z). The function E(z) in 
(48) is shown graphically in Figs. 24 and 25. 
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The successive substitution into Eq. (21) of Eqs. (20), (1), (17), and (19) gives after 
simplification 


€ 1 
ee ee ea ee 22 
iia eae aD 


a= atin. = 


This is the parabola 

y= ce (1 + Pies (23) 
Since a change in the value of c, effects only a change of scale in the ¢-plane, ¢ may 
be taken without loss of generality as 

c= eX(1 + 3, (24) 

and this parabola becomes the one considered in Eq. (2). Setting this value of ¢ in 
(19) yields 

b= 4(1 + ee. (25) 
Hence, by Eqs. (7) and (14), 


aH (lbw /i+9l+3), asm i/l—e), g=f/(ite. (26) 


Formulae (15) and (16) are valid for ¢;>0, o2>0, i.e., for e<1 and thus include those 
profiles whose thicknesses is less than about 4/5 of their lengths. 
It may also be noted that in terms of the variable y, the slope, @(y), and the 


curvature, d0(y)/ds, for the symmetrical profile may be written as 
4 tan y(tan? 1 + 2€)/e? 
6(y) = y — arctan eran ; (27) 
4 tan? y — [tan? y + (1 + 2e)/e?]? 


dé 
2 7 — [(1 + 3e)(1 — ee? cost y + 6 cos? y — 3e2/(1 + 6)*]. (28) 


CORRECTION AND SUPPLEMENT TO OUR PAPER 
THE CYLINDRICAL ANTENNA: CURRENT AND IMPEDANCE* 
QUARTERLY OF APPLIED MATHEMATICS 3, 302-335 (1946) 
’ By RONOLD KING anp DAVID MIDDLETON (Harvard University) 


Equation (58) should be written as follows: 
ee { %x:(0) | = | ¥1(0) | /sin Bh; Bhs x/2 
| Wii(h — /4)| =| yilh— 2/4) |; Bh 2 w/2. 


Two lines before this equation | ¥1(0) | /sin Bh should be written instead of | y(0)| . 

These changes involve no alternations in the figures. However, the function 
| (0) | plotted in Fig. 11 to the left of Bh=7/2 is not the parameter of expansion y 
defined by (58) as modified above and as indicated in the caption. The parameter of 
expansion y as defined in (58) is plotted in Fig. 11a where the part to the right of 
Bh=7/2 is the same as in Fig. 11, the part to the left of 8h = 7/2 is obtained from the 
curves in Fig. 11 by dividing by sin Bh. 


(58) 


* Received Jan. 25, 1946. 


200 NOTES 


For small values of Bh a convenient approximate formula is 
Vxi(0) =Q2—2—j8h; Bh<0.5 
so that 
| ¥xi(0) | = V(Q— 2)? + BH = Q— 2 + 4Bth?/(Q — 2). 


Bh 


Fic, 11a. The expansion parameter y as defined in the corrected equation (58). 


3.0 3.5 40. 


The following minor errors and misprints have been called to our attention: 

page 312, Eq. (43) change W to y, 

page 319, Eqs. (59) and 62), change Vx,(z) to y; line following Eq. (61), delete 
the following: y(z) =0 and 

page 320, Eqs. (69) and (70), change b to Q; Eq. (76), insert 1/(m—1)! after the 
first equality sign, 

page 323, Eq. (77b), change 4 to y, 

page 324, Eq. (79), insert y after R., 

page 329, Eq. (19), third line, change (Rea+u2) to (Re,—ue), 

page 330, Eqs. (23) and (27), page 335 Eqs. (45) and (46), and in the integral 
preceding Eq. (45), change Re, to uz, Ri, to %, throughout, 

page 330, Eq. (24) replace by: w=(h+2); u1=(h—z), 

page 332, Eq. (43), add superscript bar over first three symbols Ci. 


